CONTRACT  REPORT  ARBRL-CR-00528 


AD 


AERODYNAMIC  HEATING  COMPUTATIONS  FOR 
PROJECTILES  -  VOL.  II;  SWEPT  WING 
CALCULATIONS  USING  THE  PLANAR 
VERSION  OF  THE  ABRES  SHAPE 
CHANGE  CODE  (PLNRASCC) 

Prepared  by 

Acurex  Corporation,  Aerotherm  Division 
555  Clyde  Avenue,  P.  0.  Box  7555 
Mountain  View,  California  94039 


June  1984 


US  AflMY  ARMAMENT  RESEARCH  ANR  REVEEOPMENT  CENTER 

BALLISTIC  RESEARCH  LABORATORY 

ABERDEEN  PROVING  GROUND,  MARYLAND 


Approved  for  public  release;  distribution  unlimited. 


19970908  139 


3DTIG  QUALITY  mSBmiM)  ^ 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.  S.  Department  of  Commerce,  Springfield,  Virginia 

22161. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other 
authorized  documents. 


The  use  of  trade  names  or  ma>\'’A.faaturers '  names  in  this  report 
joes  not  constitute  indorsement  of  <x'.y  ootmsroidl  prcduat. 


UNCLASSIFIED _ 

security  classification  of  this  page  (Whon  Data  Entarmd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

CONTRACT  REPORT  ARBRL-CR-00528 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  title  (and  Subtitta) 

AERODYNAMIC  HEATING  COMPUTATIONS  FOR  PROJECTILES 
-  VOLUME  II:  SWEPT  WING  CALCULATIONS  USING  THE 
PLANAR  VERSION  OF  THE  ABRES  SHAPE  CHANGE  CODE 
(PLNRASCC) 

5.  type  of  REPORT  A  PERIOD  COVERED 

Final 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfa) 

Roger  C.  Strawn  and  William  S.  Kobayashi 

8.  CONTRACT  OR  GRANT  NUMBERCa; 

DAAK11-81-C-0064 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Acurex  Corporation,  Aerotherm  Division 

555  Clyde  Avenue,  P.O.  Box  7555 

Mountain  View.  California  94039  _  .. 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  a  WORK  UNIT  NUMBERS 

RDT&E  1L162618AH80 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  AMCCOM,  ARDC 

Ballistic  Research  Laboratory,  ATTN:  DRSMC-BLA-S{A) 
Aberdeen  Proving  Ground,  MD  21005 

12.  REPORT  DATE 

June  1984 

13.  NUMBER  OF  PAGES 

_ aa _ 

14.  MONITORING  AGENCV  NAME  &  AODRESS('lf  difiarant  from  Controlling  Offica) 

15.  SECURITY  CLASS,  (of  tht,  report; 

Unclassified  . 

15a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  Co/ tflla  Report; 

Approved  for  public  release,  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  tha  abatract  antarad  in  Block  20,  if  diftarant  from  Raport) 

18.  SUPPLEMENTARY  NOTES 

This  work  was  performed  under  the  direction  of  the  Aerodynamics  Research  Branch, 
Launch  and  Flight  Division,  DRSMC-BLL  (A),  Dr.  Walter  B.  Sturek,  Contracting 
Officer's  Technical  Representative. 

19.  KEY  WORDS  (Continua  on  ravaram  mida  if  nacaaaary  and  idantlfy  by  block  numbar) 


Aerodynamic  Heating 
Unsteady  Heat  Conduction 
Swept  Fin  Heating 
Numerical  Computations 

[20!  ABSTRACT  CCanttaum  mt  nvmrmm  aUt  U  imct— aptt  IdeaUtr  by  btoek  numbmr) 

This  report  documents  modifications  and  additions  incorporated  into  the 
ABRES  Shape  Change  Code  (ASCC80)  to  create  a  planar  two-dimensional  version 
of  the  axisymmetric  computer  code.  This  planar  code  predicts  convective  heat 
transfer  and  in-depth  conduction  for  swept  wings  in  supersonic  flow.  The 
report  contains  test  cases  and  a  detailed  user's  guide  which  describes  the 
input  data  required  to  run  the  code. 


DD  /jSTts  M73  EDITION  OF  f  NOV  65  |S  OBSOLETE 


TABLE  OF  CONTENTS 


Page 

LIST  OF  ILLUSTRATIONS .  5 

INTRODUCTION .  1 

TECHNICAL  DISCUSSION .  9 

2.1  Viscous  and  Inviscid  Flowfield  Modifications .  9 

2.1.1  Results . 

2.1.2  Discussion  and  Conclusions  on  Flowfield  Modifications....  21 

2.2  Planar  In-Depth  Heat  Conduction  Modifications .  22 

2.2.1  Implicit  Grid  Modifications .  22 

2.2.2  Explicit  Grid  Modifications .  23 

2.2.3  Validation  of  Planar  Conduction  Modifications .  24 

INPUT  AND  OUTPUT .  27 

3.1  Input  Modifications  to  the  BRLASCC  Computer  Code .  27 

3.2  Sample  Problems .  30 

REFERENCES . 84 

NOMENCLATURE . 85 

DISTRIBUTION  LIST . 87 


3 


LIST  OF  ILLUSTRATIONS 

Figu re  ^^9^ 

2-1  Sketch  of  Swept  Wing  Geometry .  10 

2-2  Heat  Transfer  Predictions  of  the  Experimental  Data  of 

Stainback  =  0.25  inch,  Pq  =  428  psig,  Tq  =  460°F .  12 

2-3  Heat  Transfer  Predictions  of  the  Experimental  Data  of 

Stainback  Rj  =  0.25  inch,  Pq  =  223  psig,  Tq  =  441°F .  12 

2-4  Heat  Transfer  Predictions  of  the  Experimental  Data  of 

Stainback  Rj  =  0.25  inch,  Pq  =  109  psig,  Tq  =  432°F .  13 

2-5  Heat  Transfer  Predictions  of  the  Experimental  Data  of 

Stainback  R^  =  0.25  inch,  Pq  =  65  psig,  Tq  =  437°F .  13 

2-6  Heat  Transfer  Predictions  of  the  Experimental  Data  of 
Murray  and  Stallings  R^-  =  0.125  inch,  =  3.71, 

Re  =  9.85  X  10^  per  meter,  A  =  60° .  15 

2-7  Heat  Transfer  Predictions  of  the  Experimental  Data  of 
Murray  and  Stallings  R^*  =  0.125  inch,  M„  =  3.71, 

Re  =  19.7  X  10®  per  meter,  A  =  60° . . .  15 

2-8  Heat  Transfer  Predictions  of  the  Experimental  Data  of 
Murray  and  Stallings  R-j  =  0.125  inch,  =  4.44, 

Re  =  9.85  X  10®  per  meter,  A  =  60° .  16 

2-9  Heat  Transfer  Predictions  of  the  Experimental  Data  of 
Murray  and  Stallings  Rj  =  0.125  inch,  M^  =  4.44, 

Re  =  19.7  X  10®  per  meter,  A  =  60° .  16 

2-10  Heat  Transfer  Predictions  of  the  Experimental  Data  of  Hunt, 
et  al,  R^  =  0.5  inch,  M„  =  7.81,  Re  =  0.92  x  10®  (based  on 

leading  edge  diameter),  A  =  60° . . .  18 

2-11  Heat  Transfer  Predictions  of  the  Experimental  Data  of  Hunt, 
et  al,  R^-  =  0.5  inch,  =  7.94,  Re  =  2.6  x  10®  (based  on 

leading  edge  diameter),  A  =  60° .  18 

2-12  Heat  Transfer  Predictions  of  the  Experimental  Data  of  Hunt, 
et  al ,  R^  =  0.5  inch,  =  7.98,  Re  =  9.3  x  10®  (based  on 

leading  edge  diameter),  A=  60° .  19 

2-13  Infinite  Rectangle  Configuration .  25 

2-14  Comparison  Between  PLNRASCC  Conduction  Solution  and  Heisler 

Chart  Solution  for  an  Infinite  Rectangle .  26 

5 


SECTION  1 


INTRODUCTION 

This  report  documents  modifications  and  additions  incorporated  into  the 
ABRES  Shape  Change  Code  (ASCC80)1  to  create  a  planar  two-dimensional  version 
of  the  axisymmetric  computer  code.  This  planar  code  predicts  convective  heat 
transfer  and  in-depth  conduction  for  swept  wings  in  supersonic  flow.  These 
planar  modifications  to  ASCC80  were  developed  under  the  Aerodynamic  Heating 
Computations  for  Projectiles  program.  The  overall  objectives  for  the  program 
were  threefold: 

1.  Modify  the  in-depth  heat  conduction  package  to  improve  ASCC's 
capabilities  to  handle  slender  mul timaterial  configurations 

2.  Extend  the  developments  of  planar  ASCC  modifications  to  predict 
heating  of  swept  fin  configurations  to  include:  (a)  turbulent  flow 
on  swept  wings;  (b)  2-D  shock  shape;  and  (c)  improved  in-depth  heat 
conduction  routines 

3.  Develop  an  interactive  computational  grid  develping  routine  to 
simplify  the  procedure  for  specifying  body  configurations  and 
developing  computational  grids  for  ASCC 

The  modifications  made  to  ASCC80  covering  the  second  objective  are 
documented  in  Volume  II  of  this  report.  Volumes  I  and  III  of  this  report 
document  the  work  related  to  Objectives  1  and  3,  respectively.  In  this 
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document,  the  updated  ASC  code  is  referred  to  as  PLNRASCC,  and  the  updated 
code  associated  with  Objective  1  is  referred  to  as  BRLASCC. 

Technical  discussion  of  the  PLNRASCC  modifications  is  presented  in 
Section  2.  Section  3  is  devoted  to  a  discussion  of  input  and  output. 
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SECTION  2 


TECHNICAL  DISCUSSION 

2.1  VISCOUS  AND  INVISCID  FLOWFIELD  MODIFICATIONS 

Viscous  and  inviscid  flow  models  for  the  ASCC80  version  of  the  ABRES 
Shape  Change  Code^  have  been  modified  to  solve  a  planar  two-dimensional  model 
for  flow  over  swept  wings.  The  theoretical  basis  for  these  swept  wing 
modifications  is  described  in  a  report  by  Suchsland.^  Implementation  into  the 
PLNRASCC  computer  program  is  taken  directly  from  that  report.  In  addition, 
two  new  capabilities  have  been  added  to  the  fluid  flow  modeling  in  the 
computer  code.  The  first  is  the  ability  to  model  two-dimensional  shock  shapes 
and  include  these  effects  in  the  boundary  layer  calculation.  Suchsland's 
version  of  the  code  could  only  model  shock  shapes  for  axi symmetric 
configurations.  The  second  added  capability  is  the  use  of  curve  fitted 
pressure  correlations  for  planar  geometries.  The  earlier  version  of  the 
inviscid  flow  model  was  restricted  to  axisymmetric  body  shapes. 

The  planar  shock  shape  modifications  to  the  code  are  a  series  of  curve 
fits  to  computed  results  from  a  two-dimensional  version  of  the  RAZZ I 
computer  code.  The  RAZZIB  code  is  a  general  inviscid  flow  solver  for 
highspeed  flight  configurations.  Pressure  distributions  were  computed  for  a 
series  of  cylinder-wedge  two-dimensional  wings.  Flow  conditions  for  these 
test  cases  spanned  a  range  of  Mach  numbers  (1^,  =  1.75,  2.0,  3.0,  4.0,  and 
6.0),  and  also  a  range  of  aft  wedge  angles  (0y,  =  0.0°,  2.5°,  and  5.0°).  The 
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X: 


Figure  2-1.  Sketch  of  swept  wing  geometry 


computed  results  for  shock  angle,  0y^,  versus  X/R.j  were  fitted  with  a  least 
squares  sixth  degree  polynomial  for  each  combination  of  wedge  angle  and  Mach 
number  and  implemented  into  the  BRLASC  code.  Internally,  these  curve  fits 
are  converted  into  tabular  form.  Linear  interpolation  in  both  wedge  angle  and 
Mach  number  is  used  to  create  the  shock  shape  table  for  each  new  swept  wing 
geometry.  Note  that  in  the  case  of  a  swept  wing,  the  required  shock  shape  is 
computed  using  the  component  of  freestream  Mach  number  that  is  normal  to  the 
leading  edge  of  the  wing. 

The  planar  inviscid  pressure  correlations  are  also  implemented  in  the 
form  of  curve  fits  to  inviscid  flow  calculations.  Results  from  the  RAZZIB 
code  were  obtained  for  the  same  ranges  of  Mach  numbers  and  wedge  angles 
described  above.  Curve  fits  for  P/Pq  versus  X/R-j  were  produced  using  a 
combination  of  sixth  degree  polynomial  and  exponential  least  squares 
functions.  Linear  interpolation  is  used  in  both  wedge  angle  and  Mach  number 
in  order  to  produce  a  pressure  ratio  for  each  boundary  layer  integration 
point.  As  was  noted  above  for  the  shock  shape  predictions,  it  is  the  swept 
wing's  normal  component  of  Mach  number  that  determines  where  to  interpolate  in 
the  pressure  curve  fits. 
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The  new  planar  pressure  correlations  are  only  used  on  the  aft  wedge 
portion  of  the  wing.  The  existing  ASCC80  pressure  correlations  use  modified 
Newtonian  theory  to  predict  the  pressure  distribution  on  the  nosetip  of  an 
axisymmetric  configuration.^  Newtonian  theory  applies  to  planar  geometries  as 
well  as  axisymmetric  ones.  Thus  the  original  pressure  correlations  were  not 
altered  in  the  nosetip  region. 

Planar  input  modifications  to  the  BRLASC  code  are  described  in 
Section  3.  This  section  should  be  used  in  addition  to  the  user's  manual  in 
Volume  I  of  this  report  in  order  to  run  the  new  code  for  planar  swept  wings. 
These  input  modifications  enable  PLNRASCC  to  be  considerably  more  versatile 
and  easier  to  use  than  Suchsland's  original  version  of  the  code. 

2.1.1  Resul ts 

This  section  gives  convective  heat  transfer  results  from  PLNRASCC  and 
compares  them  with  experimental  data  from  a  number  of  different  swept  wing 
configurations.  These  experiments  were  chosen  to  test  the  code  for  its 
ability  to  compute  heat  transfer  in  both  laminar  and  turbulent  boundary 
layers. 

Figures  2-2  through  2-5  show  predictions  of  the  swept  wing  data  of 
Stainback.5  This  experiment  consisted  of  a  60°  swept  delta  wing  in  a 
supersonic  flow  at  several  different  tunnel  stagnation  pressures.  The 
freestream  Mach  number  for  these  data  is  4.95.  Laminar  flow  conditions  exist 
for  the  entire  run  length  of  the  wing.  The  "new  prediction"  in  Figures  2-2 
through  2-5  was  produced  using  PLNRASCC.  The  "old  prediction"  is  taken  from 
Suchsland.2  The  differences  between  these  two  codes  lie  in  the  formulation  of 
shock  shape  and  inviscid  pressure  results  for  planar  geometries.  In 
Suchsland's  prediction,  he  assumes  a  normal  shock  at  the  leading  edge  of  the 
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o  STAINBACK  DATA 
-  OLD  PREDICTION 

- new  prediction 


0-00  1.00  2.00  3.00  4.00 

NORMALIZED  RUN  LENGTH  (s/r) 


5.00 


Figure  2-2.  Heat  transfer  predictions  of  the  experimental  data  of  StainbackS 


Ri  =  0.25  inch,  Pq  =  428  psig,  Tq  =  460°F 


O  STAINBACK  DATA 

-  OLD  PREDICTION 

-  NEW  PREDICTION 


j-  0.00  1.00  2.00  3.00  4.00  5.00 

NORMALIZED  RUN  LENGTH  (s/r) 

Figure  2-3,  Heat  transfer  predictions  of  the  experimental  data  of  Stainback^ 
Ri  =  0.25  inch,  Pq  =  223  psig.  To  =  441°F 


HEAT  FLUX  C BTU/FT*¥2-SEC )  ^  HEAT  FLUX  ( BTU/F 

0.00  4.00  S.QQ  12.0  0.00  2.00 


NORMALIZED  RUN  LENGTH  (s/r) 

2-4.  Heat  transfer  predictions  of  the  experimental  data  of  Stainback^ 
Ri  =  0.25  inch,  Pq  =  109  psig,  Tq  =  432°F 


NORMALIZED  RUN  LENGTH  (s/r) 


Figure  2-5.  Heat  transfer  predictions  of  the  experimental  data  of  Stainback^ 
Ri  =  0.25  inch,  Pq  =  65  psig,  Tq  =  437®F 


wing.  Also,  the  pressure  ratio  (P/Po)5  is  obtained  from  the  axisymmetric 
ASCC80  correlation. 

The  results  of  Figures  2-2  through  2-5  indicate  that  PLNRASCC  and 
Suchsland's  version  of  ASCC  produce  very  similar  results  for  these  cases. 

There  are  two  reasons  for  this.  First,  the  normal  component  of  Mach  number  to 
the  wing  leading  edge  is  relatively  low,  (2.47).  This  means  that  the 
influence  of  the  shock  shape  on  the  boundary  layer  will  also  be  small.  For 
higher  Mach  numbers,  the  shock  shape  will  have  a  more  significant  effect  on 
entropy  swallowing  in  the  boundary  layer.  The  second  reason  for  the 
similarity  in  predictions  is  that,  for  zero  wedge  angle,  the  pressure 
predictions  for  the  axisymmetric  and  two-dimensional  correlations  will  be 
similar.  This  is  not  necessarily  the  case  for  higher  wedge  angles. 

Figures  2-6  through  2-9  show  predicted  and  experimental  results  for  the 
data  of  Murray  and  Stallings.^  Their  experiment  tested  both  60^  and  70°  swept 
wings  in  a  wind  tunnel  with  a  range  of  freestream  Mach  numbers  and  tunnel 
stagnation  pressures.  The  aft  wedge  angle  was  zero  for  all  of  these  cases.  A 
boundary  layer  trip  was  used  at  S  =  0.637  cm  from  the  wing  leading  edge  in 
order  to  provide  for  turbulent  flow  over  the  wing. 

For  the  70°  swept  wings  tested  by  Murray  and  Stallings,  the  normal 
component  of  Mach  number  to  the  wing  leading  edge  is  too  low  (<1.75)  to  be 
covered  by  the  pressure  and  shock  shape  correlations  described  in  Section  2.1 
of  this  report.  Therefore,  predicted  results  are  only  presented  for  cases 
with  a  60°  sweep  angle.  These  predictions  were  made  with  PLNRASCC.  The 
calculation  was  done  with  specified  transition  to  turbulent  flow  at  the 
boundary  layer  trip.  Transitional  heating  was  modeled  in  this  calculation  by 
specifying  NREYCR  =  +4  in  Input  Table  1  of  the  input  data. 
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Figure  2-6.  Heat  transfer  predictions  of  the  experimental  data 
of  Murray  and  Stall ing6  R^-  =  0.125  inch,  NU  =  3.71, 
Re  =  9.85  X  10^  per  meter,  A  =  60° 
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Figure  2-7.  Heat  transfer  predictions  of  the  experimental  data 
of  Murray  and  Stall ings^  =  0.125  inch,  M»  =  3.71, 
Re  =  19.7  X  10®  per  meter,  A  =  60° 
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Figure  2-8.  Heat  transfer  predictions  of  the  experimental  data 
of  Murray  and  Stallings^  R-j  =  0.125  inch,  H„  =  4.44, 
Re  =  9.85  X  10®  per  meter,  A  =  60° 
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Figure  2-9. 
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Heat  transfer  predictions  of  the  experimental  data 
of  Murray  and  Stallings^  R-j  =  0.125  inch,  =  4.4 
Re  =  19.7  X  10^  per  meter,  A  =  60° 


Predicted  results  cannot  be  expected  to  agree  with  the  data  near  the 
boundary  layer  trip  because  the  boundary  layer  requires  a  finite  length  after 
the  trip  until  it  becomes  fully  turbulent.  This  boundary  layer  recovery 
distance  is  clearly  seen  in  the  experimental  data  of  Figures  2-6  through  2-9. 
The  computer  model,  on  the  other  hand,  can  change  abruptly  from  laminar  to 
turbulent  flow. 

Agreement  between  the  PLNRASCC  prediction  and  the  =  3.71 
experimental  data  in  Figures  2-6  and  2-7  is  similar  to  what  is  found  by  Murray 
and  Stallings.^  They  present  predicted  results  that  were  obtained  from  the 
"strip  theory"  method  of  Van  Driest.^  For  the  higher  Mach  number  cases  in 
Figures  2-7  and  2-8,  PLNRASCC  underpredicts  the  heat  transfer  on  the  aft 
portion  of  the  wing. 

Murray  and  Stallings®  estimate  the  uncertainty  in  their  measured  heat 
transfer  coefficients  to  be  10  percent  for  h  >  306  J/m^-K,  15  percent  for 
20  J/m^-s-K  <  h  <  306  J/m^-s-K,  and  20  percent  for  h  <  20  J/m^-s-K.  The 
predicted  heat  transfer  results  from  PLNRASCC  are  close  to  being  within  these 
uncertainty  ranges.  Changes  introduced  to  the  experimental  bounJary  layer  by 
the  presence  of  the  trip  may  account  for  the  remainder  of  the  discrepancies. 

It  should  be  noted  that  in  all  cases,  the  basic  shapes  of  the  PLNRASCC  heat 
transfer  predictions  show  good  agreement  with  the  experimental  data. 

A  final  comparison  between  PLNRASCC  an<'  experimental  data  is  shown  in 
Figures  2-10  through  2-12.  These  experiments  were  conducted  by  Hunt  et  al 
They  used  temperature-sensitive  paint  on  a  60°  swept  wing  in  order  to  measure 
the  heat  transfer  rates.  High  uncertainties  are  usually  associated  with  this 
phase  change  paint  technique,  although  no  specific  values  are  provided  by  Hunt 
et  al.®  The  aft  wedge  angle  was  zero  for  all  of  their  experiments.  These 
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Figure  2-10.  Heat  transfer  predictions  of  the  experimental  data  of  Hunt 
Q  et  al.«  R,-  =  0.5  inch,  =  7.81, 

2  Re  =  0.92  X  105  (based  on  leading  edge  diameter),  A  =  60° 


Figure  2-11.  Heat  transfer  predictions  of  the  experimental  data  of  Hunt 
et  al  .8  =  0.5  inch,  K.  =  7.94, 

Re  =  2.6  X  lOS  (based  on  leading  edge  diameter),  A  =  60° 
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Figure  2-12.  Heat  transfer  predictions  of  the  experimental  data  of 
Hunt  et  al Rj  =  0.5  inch,  =  7.98, 

Re  =  9.3  X  105  (based  on  leading  edge  diameter),  A  =  60° 


experiments  had  freestream  Mach  numbers  close  to  8.0  and  a  range  of  Reynolds 
numbers. 

Results  from  the  low  Reynolds  number  case  in  Figure  2-10  are  in  fairly 
good  agreement  with  the  experimental  data.  The  largest  discrepancies  occur 
near  the  leading  edge  of  the  wing.  Laminar  flow  was  observed  experimentally 
over  the  entire  wing,  and  the  PLNRASCC result  represents  a  laminar  boundary 
layer  calculation. 

Figure  2-11  gives  heat  transfer  results  for  a  higher  Reynolds  number 
case.  Again,  the  PLNRASCC  prediction  is  in  fairly  good  agreement  with  the 
data.  Laminar  flow  was  observed  experimentally  over  the  full  length  of  the 
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wing,  so  a  laminar  boundary  layer  calculation  was  performed  with  the  PLNRASC 
code. 

Heat  transfer  results  for  the  highest  Reynolds  number  case  of  Hunt 
et  al.®  are  shown  in  Figure  2-12,  Turbulent  flow  was  observed  on  the  surface 
of  the  wing,  but  it  is  difficult  to  determine  exactly  where  boundary  layer 
transition  takes  place.  Hunt  et  al concluded  from  their  data  that  the  flow 
was  turbulent  on  the  cylindrical  nose  of  the  wing.  It  apparently  laminarized 
as  it  expanded  around  the  nose  and  onto  the  flat  plate.  This  laminarized 
boundary  layer  then  went  through  a  transition  to  turbulent  flow  at 
approximately  S  =  1.2  inches  on  the  flat  plate  portion  of  the  wing. 

The  PLNRASC  code  has  no  mechanism  for  modeling  this  sort  of  a 
laminarizing  boundary  layer.  An  attempt  was  made,  however,  to  test  the 
PLNRASC  code  against  this  experimental  data.  Two  separate  computer 
predictions  were  made.  The  first  was  a  fully  laminar  calculation.  The  second 
prediction  used  a  fixed  transition  location  at  S  =  0.2  inch  from  the  leading 
edge  of  the  wing.  Figure  2-12  indicates  that  the  turbulent  prediction  shows 
better  agreement  with  the  experimental  data,  particularly  in  the  region  away 
from  the  nose  of  the  wing.  Neither  prediction  does  well  at  the  leading  edge 
of  the  wing. 

It  should  be  noted  that  a  fully  turbulent  PLNRASCC  calculation  was  also 
attempted  for  this  case.  Abrupt  transition  to  turbulence  was  specified 
immediately  after  the  laminar  series  solution  that  starts  off  the  calculation. 
This  caused  the  boundary  layer  properties  influence  coefficient  on  Stanton 
number  to  change  abruptly.  The  effect  of  this  coefficient  is  very  important 
at  the  leading  edge  of  the  wing.  The  result  was  that  the  turbulent  heat 
transfer  rate  dropped  below  the  corresponding  laminar  heat  transfer  rate  for 
the  first  few  integration  points.  This  was  clearly  an  unrealistic  result,  so 
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the  transition  location  was  moved  farther  back  on  the  nose  until  this  problem 
was  no  longer  encountered. 

Hunt  et  al.®  also  described  a  calculation  method  that  they  used  to 
predict  the  heat  transfer  rates  in  their  experiments.  This  finite  difference 
calculation  method  included  a  spanwise  momentum  equation  for  predicting  swept 
wing  cases  with  large  crossflows.  The  PLNRASCC  heat  transfer  predictions  in 
Figures  2-10  and  2-11  agree  quite  well  with  the  predictions  that  Hunt  et  al . 
present  in  their  paper.  For  the  highest  Reynolds  number  case  shown  in 
Figure  2-12,  the  laminar  PLNRASCC  prediction  is  in  good  agreement  with  a 
laminar  prediction  that  is  given  in  the  Hunt  et  al paper, 

2.1.2  Discussion  and  Conclusions  on  Flowfield  Modifications 

For  laminar  flow  at  relatively  low  supersonic  Mach  numbers,  PLNRASCC 
does  an  adequate  job  of  predicting  heat  transfer  on  swept  wings.  This  is 
evidenced  in  the  results  of  Suchsland^  and  also  in  the  results  presented  in 
Figures  2-2  through  2-5  of  this  report.  The  limitations  of  this  method  for 
laminar  flow  are  guided  by  the  basic  assumptions  that  were  used  in  the 
formulation  of  the  swept  wing  integral  equations  of  Suchsland.^ 

The  PLNRASC  code  has  been  demonstrated  to  be  capable  of  predicting 
turbulent  boundary  layer  heat  transfer  over  swept  wings  in  supersonic  flow. 

The  predicted  results  shown  in  Figures  2-6  through  2-12  appear  to  be  quite 
reasonable  considering  the  uncertainties  associated  with  the  experimental 
data. 

One  area  that  needs  to  be  investigated  further  is  the  modeling  of 
boundary  layer  transition  on  swept  wings.  Transition  criteria  based  on 
momentum  thickness  Reynolds  number  cannot  be  easily  applied  to  PLNRASCC  for 
swept  wing  cases  with  cross-flow.  This  is  because  PLNRASCC  computes  the  flow 
in  a  two-dimensional  direction  that  is  normal  to  the  leading  edge  of  the  wing. 


21 


Reynolds  numbers  are  formed  in  the  code  by  using  a  vector  component  of 
boundary  layer  edge  velocity.  Thus,  the  boundary  layer  momentum  thickness 
Reynolds  number  that  is  calculated  by  PLNRASCC  has  a  questionable  physical 
interpretation  regarding  boundary  layer  transition.  It  is  suggested  that 
users  of  PLNRASCC  supply  a  fixed  transition  location  to  the  code  that  is 
obtained  from  a  procedure  that  is  appropriate  to  three-dimensional  boundary 
layers. 

2.2  PLANAR  IN-DEPTH  HEAT  CONDUCTION  MODIFICATIONS 

All  of  the  modifications  made  to  BRLASCC  have  been  incorporated  into 
PLNRASCC,  and  the  documentation  for  these  changes  can  be  found  in  Volume  I  of 
this  report.  However,  BRLASCC  is  an  axisymmetric  shape  change  code;  hence,  the 
heat  conduction  equation  formulations  used  in  PLNRASCC  must  be  modified  to 
correctly  model  the  planar  nature  of  the  problems  of  interest. 

2.2.1  Implicit  Grid  Modifications 

The  conduction  equation  in  the  moving  orthogonal  coordinate  system 
(implicit  grid)  under  the  axisymmetric  assumption  (d/dr  =  0)  is: 


where 

Cp  =  specific  heat 

rQ  =  body  circumferential  radius  of  curvature 
rb  =  ro  +  r«cos(0) 

r^-  =  local  streamwise  radius  of  curvature 
K  =  thermal  conductivity 
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p  =  density 

n  =  surface  normal  recession  rate,  n  =  -r 
T  =  temperature 
t  =  time 

0  =  angle  between  normal  to  local  surface  and  axis  of  symmetry 
s  =  streamwise  distance  along  body 

r  =  distance  normal  to  body  surface  at  s,  measured  from  surface 
Under  the  planar  assumption,  the  conduction  equation  takes  the  form: 


The  finite-difference  equations  in  PLNRASCC  have  been  modified  to 
reflect  the  change  in  the  differential  equation,  setting  rt  to  unity 
wherever  the  body  radius  of  curvature  appears. 

2.2.2  Explicit  Grid  Modifications 

Making  use  of  the  axisymmetric  nature  of  the  problems  of  interest  the 
conduction  equation  in  BRLASCC  utilizes  a  fixed  cylindrical  coordinate 
system  (explicit  grid)  and  is  given  by: 


PLNRASCC  however  is  not  axisymmetric,  and  the  conduction  equation  takes 
the  form  of  the  two-dimensional  Cartesian  coordinate  heat  conduction  equation 
and  is  given  by: 
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P  3T  _  3  /  3T\  ,3  /  3T\ 

3t  3x  3x /  3y  3y ) 


(4) 


The  finite-difference  equations  in  PLNRASCC  have  been  modified  to 
reflect  the  change  in  the  differential  equation. 

2.2.3  Validation  of  Planar  Conduction  Modifications 

Analytical  solutions  presented  in  graphical  form  are  available  from 
Heisler^  for  the  transient  temperature  distribution  of  multidimensional 
systems. The  solution  for  the  transient  temperature  distribution  of  an 
infinite  copper  rectangle  4  inch  x  8  inch  was  used  to  check  the  planar 
conduction  modifications  made  in  PLNRASCC.  Figure  2-13  illustrates  the 
problem. 

The  rectangle  is  at  a  uniform  initial  temperature  of  530°R.  The 
surface  temperature  is  suddenly  raised  to  800®R.  The  temperature  versus  time 
at  Point  A  was  calculated  using  the  Heisler  charts  and  compared  with  PLNRASCC. 
The  results  are  shown  in  Figure  2-14.  The  agreement  between  PLNRASCC  and  the 
Heisler  chart  solution  is  excellent.  Only  at  5  s  does  the  solution  from 
PLNRASCC  depart  significantly  from  the  Heisler  chart  solution.  This  is  due  to 
the  fact  that  the  Fourier  modulus  at  5s  is  less  than  0.2,  where  the  Heisler 
charts  become  invalid. 
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Plane  of  symmetry 


Figure  2-13.  Infinite  rectangle  configuration 
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SECTION  3 


INPUT  AND  OUTPUT 

3.1  INPUT  MODIFICATIONS  TO  THE  BRLASCC  COMPUTER  CODE 

Two  of  the  input  tables  to  the  BRLASC  code  require  modification  in 
order  to  run  the  planar  version  of  the  code  for  swept  wings.  The  conventional 
forms  of  these  tables  are  described  in  Volume  I  of  this  report.  The 
modifications  that  are  described  below  should  be  implemented  in  conjunction 
with  this  reference. 

TABLE  7:  Surface  Data 

The  modification  to  this  table  consists  of  an  additional  surface  table 
that  describes  the  geometry  of  the  swept  wing  and  also  the  specified  option 
for  computing  the  inviscid  flowfield.  Whenever  a  swept  wing  configuration  is 
to  be  computed,  the  following  subtable  must  be  included  in  the  surface  data. 
This  subtable  takes  the  following  form: 


Card  No. 

Columns 

Format 

Data 

Units 

k 

1  -  5 

15 

Enter  4  (Swept  wing  configuration 
subtable) 

— 

6-10 

15 

Option  for  pressure  calculation 
=  1  :  Use  2-D  cylinder/wedge 
correlations 

^  1  :  Do  not  use  2-D  cylinder/ 
wedge  correlations 

k+1 

1  -  10 

F10.5 

Sweep  angle  of  the  wing  measured 
from  a  normal  to  the  flow  direction 

degree 
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11  -  20  F10.5 


Aft  wedge  angle 


degree 


k+2  1-5  15  Enter  -1  if  other  surface  tables 

follow;  +1  for  last  surface  table 


Note  that  if  the  2-D  cylinder/wedge  correlations  for  pressure  are  not  desired, 
the  user  has  a  number  of  options  for  specifying  pressure  that  are  described  in 
Volume  I  of  this  report.  These  include  two  ways  of  specifying  tabular  values 
for  pressure  as  well  as  using  the  axi symmetric  pressure  correlations  that  are 
built  into  BRLASCC  and  ASCC80. 

TABLE  8:  Shock  Shape  Data 

This  table  gives  the  user  the  option  of  whether  or  not  to  use  the  2-D 
shock  shape  correlations  for  swept  wing  configurations.  The  modification  for 
input  to  this  table  is  the  addition  of  another  shock  shape  flag.  This 
modification  is  described  below. 


Card  No. 

Columns 

Format 

Data 

Uni  ts 

1 

1  -  2 

12 

Enter  08  (table  number) 

-- 

2 

1  -  5 

15 

ISHFLG  —  shock  shape  flag 

— 

1-  Shock  angle  given  as  function 
of  y  coordinate 

2-  Shock  angle  given  as  function 
of  dimensionless  y  coordinate 
(y/R  ).  where  R  is  the  nose 
radius. 

3-  Flag  to  use  the  2-D  planar  shock 
shape  correlations 


The  remainder  of  this  table 
report.  Note  that  if  TABLE 
code  will  use  a  shock  shape 
correlations  that  are  built 


is  identical  to  the  one  in  Volume  I  of  this 
8  is  not  specified  in  the  program  input  data,  the 
that  is  generated  from  the  axi symmetric 
in  to  BRLASCC  and  ASCC80. 
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SAMPLE  PROBLEM  1 


Swept  Wing  Configuration 
Boundary  Layer  Solution 
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3.2  SAMPLE  PROBLEM 


This  subsection  illustrates  two  sample  cases  for  the  PLNRASC  code.  The 

5 

first  sample  case  is  a  prediction  of  a  swept  wing  experiment  from  Stainback. 
Heat  transfer  results  from  this  case  are  presented  in  Figure  2-1. 

The  complete  input  data  for  this  case  is  presented  followed  by  selected 
portions  of  the  output  file.  The  user  should  note  that  the  boundary  layer 
calculation  proceeds  in  a  direction  that  is  normal  to  the  leading  edge  of  the 
wing.  Output  quantities  such  as  boundary  layer  edge  velocity  represent  vector 
components  of  the  total  flowfield  quantities. 

The  second  sample  problem  is  a  planar  calculation  of  the  sample  problem 
found  in  Volume  I  of  this  report.  The  axisymmetric  projectile  of  Volume  I  has 
been  made  into  a  planar  60®  swept  wing.  A  complete  listing  of  the  input  data 
is  presented  as  well  as  selected  output  data. 
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PLNRASCC  Sample  Problem  1  Input  Data 


BRL/ASCC  TWO-DIMENSIONAL  TEST  CASE 
PLANAR  SWEPT  WING 


8/10/83 


01 

0.0  0.0 

2  2  2 

1.0 

01  1.0 

02 

0.0 
1.0 
2.0 
02  3.0 
03 

50  20  1 

0.25 
04 

4.0 
4.5 
5.0 
04  5.5 

07 

4 

60.00 


2  0  1 
2.0 
4.0 

65.00  437.0 

109.0  432.0 

223.0  441.0 

428.0  460.0 


2.0  .0010 

470.0 

503.0 

537.0 

575.0 


0.0 


08 

-1 


3 

-1 


0 


4.95 


-1.00  550. 
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Sample  Problem  1  Output 


PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


4i  *  «  *  «  *  4t  *  *  *  *  «  * 


BRL/ASCC  TWO-DIMENSIONAL  TEST  CASE 
PLANAR  SWEPT  WING 


8/10 


-  GENERAL  PROGRAM  FLAGS  - 


(ENVIRONMENT  FLAG) 

(SHAPE  CHANGE  FLAG) 

(OUTPUT  PRINT  FUG) 

(TRANSITION  CRITERIA  FLAG) 
(BODY  ANGLE  DEFN.  FLAG) 

(CARBON  TRANS.  CRIT.  FLAG) 
(NOSE  SHAPE  MODIFICATION  FLAG) 


LG  =  2 

ISS  =  2 

IPRNT  =  2 

NREYCR  =  2 
IRON  =  0 

ICARB  =  1 

IMOD  =  0 


-  TIME  INCREMENT  INFORMATION  - 


INITIAL  TIME  (SEC)  0.0000 


FINAL  TIME  (SEC)  0.0000 


OUTPUT  INTERVAL  =  1.0000  SEC  FROM 

OUTPUT  INTERVAL  =  1.0000  SEC  FROM 


INITIAL  TIME  UNTIL  2.0000  SEC 
2.0000  SEC  UNTIL  FINAL  TIME 


BOUNDARY  LAYER  ONLY  SOLUTION 


COMPUTATION  TIME  STEPS  SET  EQUAL  TO  SPECIFIED  ENVIRONMENT  TIME  STEPS 


-  WIND  TUNNEL  ENVIRONMENT  - 

FREESTREAM  MACH  NO.  =  4,95 

TOTAL  TEMPERATURE 
(DEG  F) 
437.00 
432 . 00 
441 .00 
460 . 00 


TIME 
(SEC) 
0.000 
1  .000 
2.000 
3.000 


TOTAL  PRESSURE 
(PSIA) 
65.00 
109.00 
223.00 
428 . 00 
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SAMPLE  PROBLEM  2 


Swept  Wing  Configuration  With 
Transient  Heat  Conduction 


Sample  Problem  2  Input  Data 


0  0  0  0 

BRL  FLIGHT  CASE  (YUMA  TS=125  DEG-F,  T0=60  DEG-F) 

TRANSIENT  CONDUCTION  SOLUTION  —  PLNRASCC  05  JANUARY  1984 

12.5  DEG  NOSE.  7  INCH  BODY  <  BRL  SAMPLE  PLANAR  CASE>  PLNR  TEST 
01  Program  Constants  and  Time  Information 


0,0 

2.00 

4  0 

1  5  0 

1 

0.01 

0.01 

0.25 

0.25 

01 

0.25 

2.00 

02 

Environment  Table 

0.0 

1  .0 

520. 

5259 

0.2 

1  .0 

520. 

5184 

0.4 

1  .0 

520. 

5082 

0.6 

1  .0 

520. 

4941 

0.8 

1.0 

520. 

4793 

1  .0 

1  .0 

520. 

4636 

1.2 

1  .0 

520, 

4446 

1.4 

1.0 

520. 

4249 

1 .6 

1.0 

520. 

4035 

1  .8 

1.0 

520. 

3839 

02 

2.0 

1.0 

520. 

3629 

03 

Surface 

Geometry  and 

Grid  Size 

-13 

0.2  1-5  585.  75. 


0.4507 

0.0 

1 

0.4627 

0.0684 

1 

0.4975 

0.1286 

1 

0.6074 

0.1953 

1 

0.8 

.2380 

1 

1  .0 

.2823 

1 

1.249 

.3375 

1 

1.55 

.3798 

1 

1  .9 

.429 

2 

3.0 

.5836 

2 

4.0 

.7241 

2 

5.5 

0.9349 

2 

7.0 

1.1458 

2 

34 

0.4507 

0.0 

1 

0.4627 

0.0684 

1 

0.4975 

0.1286 

1 

0.6074 

0.1953 

1 

0.8 

.2380 

1 

1.0 

.2823 

1 

1.249 

.3375 

1 

1  .55 

.3798 

1 

1  .875 

.425 

1 

1  .94 

.32 

1 

1.76 

.291 

1 

1  .76 

.226 

1 

1.268 

.1571 

1 

1.268 

.0775 

1 

0.80 

.0775 

1 

0.80 

0.0 

1 
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0.4507  0.0  1 

0.80  0.0  2 

0.80  .0775  2 

1.76  .0775  2 

1.76  .291  2 

1.94  .32  2 

1.875  .425  2 

2.0  .4431  2 

4.0  .7241  2 

6.0  1.0052  2 

7.0  1.1458  2 

7.0  0.0  2 

0.80  0.0  2 

1.268  .0775  3 

1.268  .1571  3 

1.76  .2260  3 

1.76  .0775  3 

1.268  . 0775  3 

03  7 . 0  0,0 

11 
.1 
.12 

60  22 


0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1 186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.1186 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

0.055 

06  Material 

Prope  rt i es 

Table 

1  0 

1 

10.0 

0.0 

418.08 

536.0 

0.0 

400. 

.083 

1 . 794E-02 

.15 

500. 

.088 

1 .805E-02 

.15 

600. 

.093 

1 .790E-02 

.15 

700. 

.098 

1 .758E-02 

.15 

800. 

.103 

1.718E-02 

.15 

900. 

.108 

1 . 670E-02 

.15 

1000. 

.112 

1 .618E->02 

.  15 

1100. 

.117 

1 .563E-02 

.15 

-1  1200. 

.122 

1 .506E-02 

.15 

2  0 

1 

10.0 

0.0 

490. 

536.0 

0.0 
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492. 

.  1  1 

7.360E'-03 

60 

672. 

.11 

7.220E-02 

60 

1032. 

.11 

6.940E-02 

60 

-1  1392. 

.11 

6.11 0E-03 

60 

3 

0 

1 

0.0 

0.0 

488.8 

536. 

117.0 

3310. 

50. 

.128 

.1932 

540. 

.128 

7.56E-03 

5 

3250. 

.128 

7.56E-03 

5 

3320. 

.1932 

2. 10E-03 

5 

+1  9000. 

.1932 

2. 10E-03 

5 

1 

2  1  . 

0E-05 

1 

3  1  . 

0E-04 

2 

3  1.0E-06 

07  Surface 

Pressure 

Table 

4 

1 

60. 

5. 

+1 

08  Shock  Shape  Table 

yj 

09  Surface 

(  Chem i st  ry 

'  Tables 

1 

1  .0 

1  .00 

.00 

99.00 

662.210  .00 

99.556 

58-744 

1  AC'ilA- 

1  .00 

.00 

10  00 

659 . 444  . 00 

89.850 

58 . 408 

1 

1  .00 

.00 

1  .00 

656.666  .00 

89.128 

58.087 

1  AC -IK* 

1  .00 

.00 

.01 

653 . 888  . 00 

88.405 

57.765 

1  AC'^-iA^ 

1  .00 

.00 

.001 

600 . 000  . 00 

74.555 

74,555 

1  AIR 

1  .00 

.00 

.0001 

298.000  .00 

0. 

0. 

1  AIR 

1  .00 

.00 

. 00001 

200.000  .00 

-42.130 

-42.130 

1  Alt:; 

2 

1  .0 

1  .00 

.00 

.01 

2500.00  .00 

525.096 

525.096 

1  AIR 

1  .00 

.00 

.001 

600.000  .00 

74.555 

74.555 

1  AIR 

1.00 

.00 

.0001 

298.000  .00 

0. 

0. 

1  A’R 

1.00 

.00 

.00001 

200.000  .00 

-42.130 

-42.130 

1  A»^ 

3 

1.0 

1  .00 

.00 

.01 

2500.00  .00 

525.096 

525.096 

1  A  I R 

1.00 

.00 

.001 

600 . 000  . 00 

74.555 

74.555 

1  AIR 

1.00 

.00 

.0001 

298.000  .00 

0. 

0. 

1  AIR 

1  .00 

.00 

.00001 

200.000  .00 

-42.130 

-42.130 

1  AIR 

. 990+02 
.100+02 
. 100+01 
. 100-01 


-1 


-1 


Sample  Problem  2  Output 

PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


INPUT** 


BRL  FLIGHT  CASE  (YUMA  TS«125  DEG^F,  T0=60  OEG-F) 

TRANSIENT  CONDUCTION  SOLUTION  —  PLNRASCC  05  JANUARY  1984 

12.5  DEG  NOSE.  7  INCH  BODY  <  BRL  SAMPLE  PLANAR  CASE>  PLNR  TEST 


-  GENERAL  PROGRAM  FLAGS  - 


(ENVIRONMENT  FLAG) 

LG  = 

4 

(SHAPE  CHANGE  FLAG) 

ISS  » 

0 

(OUTPUT  PRINT  FLAG) 

IPRNT  - 

1 

(TRANSITION  CRITERIA  FLAG) 

NREYCR  = 

5 

(BODY  ANGLE  DEFN.  FUG) 

IRON  * 

0 

(CARBON  TRANS.  CRIT.  FLAG) 

ICARB  » 

1 

(NOSE  SHAPE  MODIFICATION  FLAG) 

IMOO  » 

0 

-  TIME  INCREMENT  INFORMATION  - 


INITIAL  TIME  (SEC)  0.0000 


FINAL  TIME  (SEC)  2.0000 


OUTPUT  INTERVAL  = 
OUTPUT  INTERVAL  « 
OUTPUT  INTERVAL  * 


0.0100  SEC  FROM 
0.2500  SEC  FROM 
0.2500  SEC  FROM 


INITIAL  TIME  UNTIL 
0.0100  SEC  UNTIL 
0.2500  SEC  UNTIL 


0.0100  SEC 
0.2500  SEC 
FINAL  TIME 


TIME  STEP  STABILITY  CRITERIA  IN  EFFECT 


MINIMUM  TIME  STEP  =  1  .000E-06  SECONDS 
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PUNAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


-  GENERAL  ENVIRONMENT  - 


TIME 

PRESSURE 

TEMPERATURE 

VELOCITY 

(SEC) 

(ATM) 

(DEG  R) 

(EPS) 

0.000 

1 .000 

520.00 

5259.00 

0.200 

1  .000 

520.00 

5184.00 

0.400 

1  .000 

520.00 

5082.00 

0.600 

1  .000 

520.00 

4941 .00 

0.800 

1  .000 

520.00 

4793.00 

1.000 

1  .000 

520.00 

4636.00 

1  .200 

1  .000 

520.00 

4446 . 00 

1.400 

1  .000 

520.00 

4249.00 

1.600 

1  .000 

520.00 

4035.00 

1.800 

1  .000 

520.00 

3839.00 

2.000 

1  .000 

520.00 

3629.00 
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PLANAR  VERSION  BRL  IMPROVED  A8RES  SHAPE  CHANGE  CODE  (PLNARASCC) 


51 


PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


- IMPLICIT  NODE  SPACING - 

NODE  THICKNESS  IN  INCHES 

NODE  NO.  !1  2  3  4  5  6  7  8  9 


BODY  PT  NO. 


1 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

2 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

3 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

4 

i  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

5 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

6 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

7 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

8 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

9 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

10 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

11 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0,0120 

0.0120 

0.0120 

12 

t  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

13 

1  0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

10  11 
- 1 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 

0.0120 
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PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


-  EXPLICIT  GRID  GEOMETRY  - 

NUMBER  OF  COLUMNS  >  60 
NUMBER  OF  ROWS  -  22 


VARIABLE  GRID  SPACING  WITH 
X  GRID  SPACING  - 


0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860 

0.11860  0.11860 

Y  GRID  SPACING  - 

0.11860 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

0.05500 

INITIAL  TEMPERATURE  OF  MODEL  -  585.0  DEG  R 
MAXIMUM  DESIRED  SURFACE  TEMPERATURE  RISE  BETWEEN  TIME  STEPS  -  75.0  DEG  R 

MINIMUM  EXPLICIT  NODAL  SPACING  USED  IN  TIME  STEP  COMPUTATION  -  0.0550  INCH 
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PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 
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PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


- MATERIAL  PROPERTIES - 


••«****  material  number  1 


-  SURFACE  ROUGHNESS 


ROUGHNESS  HEIGHT  FOR  TRANSITION 
ROUGHNESS  HEIGHT  FOR  TURBULENT  HEATING 

LAMINAR  HEATING  AUGMENTATION  FLAG 


K-LAM  -  ie.000 

K-TURB  -  0.000 

J ROUGH  -  1 


Sit! 


PARAMETERS  IN  BLOWING  CORRECTION  TO  TRANSFER  COEFFICIENTS 


LAMINAR  SHEAR  PARAMETER 
LAMINAR  HEATING  PARAMETER 
TURBULENT  SHEAR  PARAMETER 
TURBULENT  HEATING  PARAMETER 


0.5000 
0 . 5000 
0 . 3500 
0 . 3500 


-  THERMAL  PROPERTIES 


MATERIAL  DENSITY 

(RHO)  « 

418.08  ( 

;lbm/ft3) 

DATUM  TEMP  FOR  HEAT  OF 

FORMATION  (TFO)  » 

536.00  ( 

DEG  R) 

HEAT 

OF  FORMATION 

(HFO)  * 

0.00  ( 

;btu/lbm) 

TEMPERATURE 

SPECIFIC  HEAT 

CONDUCTIVITY 

SENSIBLE 

EMISSIV] 

ENTHALPY 

(DEG  R) 

(BTU/LB-OEG) 

(BTU/FT-SEC-DEG) 

(BTU/LB) 

400 . 00 

0.0830 

0.0179400 

-11.81 

0. 1500 

500.00 

0.0880 

0.0180500 

-3.26 

0.1500 

600.00 

0.0930 

0.0179000 

5.79 

0.1500 

700.00 

0.0980 

0.0175800 

15.34 

0.1500 

800.00 

0.1030 

0.0171800 

25.39 

0.1500 

900.00 

0.1080 

0.0167000 

35.94 

0.1500 

1000.00 

0.1120 

0.0161800 

46.94 

0.1500 

1100.00 

0.1170 

0.0156300 

58.39 

0.1500 

1200.00 

0.1220 

0.0150600 

70.34 

0.1500 

-  EROSION  LAW  MATERIAL  FLAG  - 

NERODE  -  0 
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PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


MATERIAL  NUMBER  2 


-  SURFACE  ROUGHNESS  - 

ROUGHNESS  HEIGHT  FOR  TRANSITION  K-LAM  -  10.000  (MIL) 

ROUGHNESS  HEIGHT  FOR  TURBULENT  HEATING  K-TURB  -  0.000  (MIL) 

LAMINAR  HEATING  AUGMENTATION  FLAG  JROUGH  «  1 


-  PARAMETERS  IN  BLOWING  CORRECTION  TO  TRANSFER  COEFFICIENTS  - 

LAMINAR  SHEAR  PARAMETER  (BLS)  »  0.5000 
LAMINAR  HEATING  PARAMETER  (BLH)  »  0.5000 
TURBULENT  SHEAR  PARAMETER  (BTS)  -  0.3500 
TURBULENT  HEATING  PARAMETER  (BTH)  *  0.3500 


-  THERMAL  PROPERTIES  - 


MATERIAL  DENSITY 

(RHO)  » 

490.00 

(LBIy(/FT3) 

DATUM  TEMP  FOR  HEAT  OF  FORMATION 

(TFO)  - 

536.00 

(DEG  R) 

HEAT  OF  FORMATION 

(HFO)  ** 

0.00 

(BTU/LBM) 

TEfctf^ERATURE  SPECIFIC  HEAT  CONDUCTIVITY  SENSIBLE 

ENTHALPY 

(DEG  R)  (BTU/LB-DEG)  (BTU/FT-SEC-DEG)  (BTU/LB) 

492.00  0.1100  0.0073600  -4.84 

672.00  0.1100  0.0722000  14.96 

1032.00  0.1100  0.0694000  54.56 

1392.00  0.1100  0.0061100  94.16 


-  EROSION  LAW  MATERIAL  FLAG  - 

NERODE  =  0 


EMISSIVITY 


0.6000 
0 . 6000 
0 . 6000 
0.6000 
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PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 

«•**•**  material  number  3  ••••♦** 


-  SURFACE  ROUGHNESS  - 


ROUGHNESS  HEIGHT  FOR  TRANSITION 
ROUGHNESS  HEIGHT  FOR  TURBULENT  HEATING 


K-LAM 

k-TURB 


LAMINAR  HEATING  AUGMENTATION  FLAG 


J ROUGH 


o.eee  (mil 
e.eee  (mil 


_  parameters  in  blowing  CORRECTION  TO  TRANSFER  COEFFICIENTS  - 

LAMINAR  SHEAR  PARAMETER  (BLS)  -  0.5000 
LAMINAR  HEATING  PARAMETER  (BLH)  =  0.5000 
TURBULENT  SHEAR  PARAMETER  (BTS)  ■  0.3500 
TURBULENT  HEATING  PARAMETER  (BTH)  -  0.3500 


-  THERMAL  PROPERTIES  - 


MATERIAL  DENSITY  (RHO) 
DATUM  TEMP  FOR  HEAT  OF  FORMATION  (TFO^ 
HEAT  OF  FORMATION  (HFO) 
LATENT  HEAT  OF  FUSION  (XLATHT) 
MELT  TEMPERATURE  (TMELT) 
TEMP  DIFFERENCE  MELT  OCCURS  (DTMELT) 
SPECIFIC  HEAT  OF  SOLID  (CPSOLD) 
SPECIFIC  HEAT  OF  LIQUID  (CPLIQO) 


488.80  (LBM/FT3) 

536.00  (DEG  R) 

0.00  (BTU/LBM) 

117.00  (BTU/LBM) 

3310.00  (DEG  R) 

50.00  (R  DEG) 

1.280E-01  (BTU/LBM-DEG  R) 
1.932E-01  (BTU/LBM-DEG  R) 


TEMPERATURE 

(DEG  R) 
540.00 
3250.00 
3285.00 
3335.00 
3336.00 
9000.00 


SPECIFIC  HEAT 

CONDUCTIVITY 

SENSIBLE 

ENTHALPY 

(BTU/L8-DEG) 

(BTU/FT-SEC-DEG) 

(BTU/LB) 

0.1280 

0.0075600 

0.51 

0.1280 

0.0075600 

347 . 39 

0.1280 

0.0048300 

351 .87 

4.5520 

0.0021000 

468 . 87 

0.1932 

0.0021000 

471.24 

0.1932 

0.0021000 

1 565 . 53 

EMISSIVITY 


0.5000 
0 . 5000 
0.5000 
0.5000 
0 . 5000 
0 . 5000 


-  EROSION  LAW  MATERIAL  FLAG  - 

NERODE  =  0 
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PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


MAT1 

CONTACT 

MAT2 

RESISTANCES 

RESISTANCE 

(^♦♦2-S~DegR/BTU) 

1 

2 

1 .00000E-05 

1 

3 

1 .00000E-04 

2 

3 

1 .00000E-06 

PLANAR  SWEPT  WING  CASE 

SWEPT  WING  ANGLE  «  60.00000  (DEG) 

WEDGE  ANGLE  - 

5.00000  (DEG) 

2-D  CYLINDER-WEDGE  PRESSURE  CORRELATIONS  ARE  USED 

MACH  NUMBER  NORMAL  TO  LEADING  EDGE  >  2.34182 
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PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


DIMENSIONLESS 

Y-COORDINATE 

Y/RN 


0.000 
0.100 
0.200 
0.300 
0.400 
0.500 
0.600 
0.700 
0.800 
0.900 
1  .000 

1.100 
1.200 

1.300 

1 .400 

1.500 

1.600 

1.700 

1.800 

1.900 
2.000 

2.100 
2.200 

2.300 

2.400 

2.500 

2.600 

2.700 

2.800 

2.900 
3.000 

3.100 

3.200 

3.300 

3.400 

3.500 

3.600 

3.700 

3.800 

3.900 
4.000 

4.100 

4.200 

4.300 


SHOCK  SHAPE 
SHOCK  ANGLE 


(DEGREES) 

90 . 000 
88.772 
87 . 276 
85.819 
84.400 
83.018 
81 .673 
80.363 
79.090 
77.850 
76 . 645 
75.474 
74.335 
73.228 
72.153 
71 .109 
70.096 

69.112 
68.157 
67 . 230 
66.332 
65.460 
64.615 
63.797 
63.003 
62 , 235 
61 .491 
60 . 770 
60 . 073 
59.397 
58.744 

58.113 
57 , 502 
56.911 
56 . 340 
55.788 
55.255 
54 . 739 
54.242 
53.761 
53.297 
52.849 
52.416 
51 .999 


DIMENSIONLESS 

X-COORDINATE 

X/RN 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


4.400 

51 .596 

0.000 

4.500 

51.207 

0.000 

4.600 

50.832 

0.000 

4.700 

50.470 

0.000 

4.800 

50.120 

0.000 

4.900 

49.783 

0.000 

5.000 

49.458 

0.000 

5.100 

49.144 

0.000 

5.200 

48.841 

0.000 

5.300 

48.549 

0.000 

5.400 

48.267 

0.000 

5.500 

47.995 

0.000 

5.600 

47.732 

0.000 

5.700 

47 . 478 

0.000 

5.800 

47.233 

0.000 

5.900 

46.996 

0.000 

6.000 

46.767 

0.000 

6.100 

46.546 

0.000 

6.200 

46.332 

0.000 

6.300 

46.125 

0.000 

6.400 

45.925 

0.000 

6.500 

45.731 

0.000 

6.600 

45.543 

0.000 

6.700 

45.361 

0.000 

6.800 

45.184 

0.000 

6.900 

45.012 

0.000 

7.000 

44.846 

0.000 

7.100 

44.684 

0.000 

7.200 

44.527 

0.000 

7.300 

44.373 

0.000 

7.400 

44.224 

0.000 

7.500 

44.079 

0.000 

7.600 

43.937 

0.000 

7.700 

43.798 

0.000 

7.800 

43.663 

0.000 

7.900 

43.531 

0.000 

8.000 

43.401 

0.000 

8.100 

43.274 

0.000 

8.230 

43.113 

0.000 

8.399 

42.909 

0.000 

8.619 

42.653 

0.000 

8.904 

42.332 

0.000 

9.276 

41 .933 

0.000 

9.758 

41 .438 

0.000 

10.386 

40,827 

0.000 

11.201 

40.080 

0.000 

12.262 

39.198 

0.000 

13.641 

38.263 

0.000 

15.433 

37.455 

0.000 

17.763 

37.112 

0.000 

20.791 

36.538 

0.000 

24.729 

35.711 

0.000 

29.847 

34.637 

0.000 

36.501 

33.240 

0.000 

45.152 

31 .424 

0.000 

56.397 

30.406 

0.000 

71 .017 

30.406 

0.000 
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PLANAR  VERSION  BRL  IMPROVED  ABRES  SHAPE  CHANGE  CODE  (PLNARASCC) 


- SURFACE  EQUILIBRIUM  DATA - 


MAT  =  1 

CMH  -  1 . 00000 


MASS  TRANSFER  COEF.  -  0.0000  LBM/FT* 


TEMP 

BPRIM 

HCH 

360.0000 

0.0000 

-15.2280 

536.4000 

0.0001 

0.0362 

1080.0000 

0.0010 

56.1020 

1176.9984 

0.0100 

67 . 5933 

1181.9988 

1 . 0000 

68.1909 

1186.9992 

10.0000 

68 . 7884 

1191 .9780 

99 . 0000 

69.3834 

MAT  -  2 

CMH  -  1 . 00000 


MASS  TRANSFER  COEF.  -  0.0000  LBM/FT** 


TEMP 

BPRIM 

HCH 

360.0000 

0 . 0000 

-19.3600 

536 . 4000 

0.0001 

0.0440 

1080.0000 

0.0010 

59.8400 

4500 . 0000 

0.0100 

436.0400 

MAT 

-  3 

CMH 

-  1.00000 

MASS  TRANSFER 

COEF.  -  0.0000 

LBM/FT 

TEMP 

BPRIM 

HCH 

360 . 0000 

0 . 0000 

-22 . 5280 

536.4000 

0.0001 

0.0512 

1080.0000 

0.0010 

69 . 6320 

4500 . 0000 

0.0100 

696.1294 

■SEC 

PRESSURE  » 

1.0000  ATM 

TSEN 

TCHEM 

SPECIE 

-75.8340 

75.8346 

AIR 

0.0000 

0.0000 

AIR 

134.1990 

-134.2771 

AIR 

103.9770 

-159.4928 

AC41 

104.5566 

-196.7961 

AC41 

105.1344 

-525.1900 

AC41 

105.7392 

-3778.4279 

AC41 

-SEC 

PRESSURE  » 

1.0000  ATM 

TSEN 

TCHEM 

SPECIE 

-75.8340 

75.8346 

AIR 

0 . 0000 

0 . 0000 

AIR 

134.1990 

-134.2734 

AIR 

945.1728 

-950.2641 

AIR 

SEC 

PRESSURE  = 

1 . 0000  ATM 

TSEN 

TCHEM 

SPECIE 

-75.8340 

75 . 8345 

AIR 

0 . 0000 

0 . 0000 

AIR 

134.1990 

-134.2636 

AIR 

945.1728 

-947.6632 

AIR 
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NOMENCLATURE 


H  convective  heat  transfer  coefficient  in  Figures  2-5  through  2-11 

A  swept  wing  angle  measured  from  a  normal  to  the  flow  direction 

Mco  freestream  Mach  number 

P  static  pressure 

Pq  stagnation  pressure 

Re  Reynolds  number 

R-j  nose  radius  on  cylinder/wedge  swept  wing 

S  distance  along  surface  measured  normal  to  the  wing  leading  edge 

T  static  temperature 

Tq  stagnation  temperature 

X  axial  distance  measured  normal  to  the  wing  leading  edge 

0y^  aft  wedge  angle  on  cylinder/wedge  swept  wing 
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